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Abstract 

We consider the best today available observations 
of the Sun free of turbulent Earth atmospheric effects, 
taken with the Solar Optical Telescope (SOT) onboard 
the Hinode spacecraft. Both the instrumental smear¬ 
ing and the observed stray light are analyzed in order 
to improve the resolution. The Point Spread Function 
(PSF) corresponding to the blue continuum Broadband 
Filter Imager (BFI) near 450 nm is deduced by analyz¬ 
ing i/ the limb of the Sun and ii/ images taken during 
the transit of the planet Venus in 2012. A combination 
of Gaussian and Lorentzian functions is selected to con¬ 
struct a PSF in order to remove both smearing due to 
the instrumental diffraction effects (PSF core) and the 
large-angle stray light due to the spiders and central ob¬ 
scuration (wings of the PSF) that are responsible for the 
parasitic stray light. A Max-likelihood deconvolution 
procedure based on an optimum number of iterations 
is discussed. It is applied to several solar held images, 
including the granulation near the limb. The normal 
non-magnetic granulation is compared to the abnormal 
granulation which we call magnetic. A new feature ap¬ 
pearing for the hrst time at the extreme- limb of the 
disk (the last 100 km) is discussed in the context of the 
dehnition of the solar edge and of the solar diameter. 
A single sunspot is considered in order to illustrate how 
effectively the restoration works on the sunspot core. A 
set of 125 consecutive deconvolved images is assembled 
in a 45 min long movie illustrating the complexity of 
the dynamical behavior inside and around the sunspot. 
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1 Introduction 

A new era in solar surface imaging techniques at high 
resolution has emerged in recent years thanks to i/ 
the introduction of adaptive optics and computer aided 
speckle imaging methods at ground-based solar obser¬ 
vatories with 1 to 1.6 m diameter aperture telescopes; 
ii/ the availability of a wealth of diffraction limited im¬ 
ages from the Solar Optical Telescope (SOT) of the 
Hinode mission, (Tsuneta et ah, 2008; Suematsu et ah, 
2010) free of turbulent Earth atmospheric smearing ef¬ 
fects with perfectly reproducible properties. It is then 
useful to carefully analyze ah instrumental effects influ¬ 
encing these space-borne images in order to extract the 
hnest details of the solar surface up to its extreme limb, 
in particular for the knowledge of dynamical effects pro¬ 
duced by the magnetic helds continually emerging to 
the surface. This includes the case of very large con¬ 
centration of magnetic flux as is the case in sunspots, 
although only the case of the less complicated single, 
unipolar, well developed sunspot near the center of the 
Sun will be considered here. We leave for a another pa¬ 
per the main application we want to consider of sunspot 
evolution, without addressing the question of the origin 
of the sunspot phenomenon in deep layers of the Sun. 

A qualitative improvement in very hne details of the 
dynamical layer immediately above the photosphere, as 
observed using an H hne emissions of Call hlter, has 
already been the subject of several Hinode SOT pa¬ 
pers devoted to the study of spicules, their roots and 
their extension towards the corona (see for ex. Tavabi, 
Koutchmy and Adjabshirzadeh, 2011 & 2014). We will 
not apply the rather qualitative method described in 
Tavabi, Koutchmy and Adjabshirzadeh (2013), in the 
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case of photospheric surface features. Instead, more 
quantitative methods should be used considering phe¬ 
nomena dominated by radiative transfer aspects, where 
the optical thickness effects are important, which was 
not the case for spicules observed off-disk with a broad 
band H Call filter (Tavabi, Koutchmy and Adjab- 
shirzadeh, 2011). Telescopes cannot work in practice 
at their ultimate diffraction-limited quality due to op¬ 
tical aberration(s) arising from imperfections, thermal 
changes, dirt and maladjustments from optimum pa¬ 
rameters after working for a period of time, including 
the defocusing. In the case of the Hinode SOT we deal 
with very small effects. 

Scattered light and stray light are mathematically 
described using the extended wings of the point spread 
function (PSF), which causes an image of the point 
source of light to be spread (Johnson, 1972). This effect 
can be enhanced when long averaging exposure times 
are used in ground based observations. In ideal condi¬ 
tions of a non-obstructed pupil, the PSF of a diffrac¬ 
tion limited telescope is described by the Airy pattern. 
Usually optical telescopes work in situations that are 
far from these ideal conditions. Recall that in the case 
of a ground based instrument the earth’s atmosphere 
will limit observations. The access to seeing-free data 
should be of great benefit with space-born instruments. 
For space-born telescopes, the PSF was tentatively de¬ 
termined in laboratory conditions before the flight (Sue- 
matsu et ah, 2010). However, after the launch some 
in-flight adjustments are performed and because the 
whole instrument works in a different situation in space, 
it is better to try to measure the PSF ”in situ”, us¬ 
ing indirect methods such as the Mercury transit, solar 
eclipses, or using the limb of the Sun (see Mathew et 
ah, 2009; Wedemeyer-Bohm, 2008; Wedemeyer-Bohm 
and Rouppe van der Voort, 2009), for earlier deductions 
of the PSF of the Hinode SOT. Note that its optical 
telescope assembly (OTA) has a 50 cm outer diame¬ 
ter with three spiders (4 cm wide) holding the 17.2 cm 
central mirror, producing an obstruction in the path of 
light (Tsuneta et ah, 2008; Suematsu et ah, 2008).The 
diffraction-limited PSF of a circular aperture of 50 cm 
clear aperture is described by the Airy function where 
the wavelength is a parameter. For more complicated 
apertures the PSF can be deduced from the Fourier 
transform of the pupil shape when it is known. Note 
that the central obscuration (in linear scale 0.344) is 
not exactly a part of the pupil because it is well before 
the main mirror and it is illuminated by the full Sun 
’’from the back”, producing some amount of spurious 
stray light towards the primary mirror that is further 
reduced by the Gregorian system used in the telescope 
(see Suematsu et ah, 2008). The theoretical diffraction- 
limited PSF with three extended spiders was discussed 


in Wedemeyer-Bohm and Rouppe van der Voort (2009) 
in terms of the modulation transfer function (MTF), 
which is the Fourier transform of the PSF. 

It is not easy to apply the concept to image restora¬ 
tion using the ’’theoretical” MTF because the effect 
of the noise dominates at high spatial frequencies and 
should be taken into account. The Wiener 2D filter 
(see Koutchmy and Koutchmy, 1989 and Bracewell, 
1976) is a statistically justified solution but it does not 
take into account the local conditions with a position- 
dependent noise. Wedemeyer-Bohm (2008) utilized an 
eclipse and the Mercury transit to evaluate the PSF 
of the SOT by qualitatively introducing a convolution 
of an ’’ideal” PSF part (Airy function) and a ’’non¬ 
ideal” PSF part (overall, described by a Voigt func¬ 
tion). Note that the edge of the Moon is far from be¬ 
ing smooth, with the profile produced by mountains 
and craters to be reconstructed for the time of obser¬ 
vation because of the lunar libration and this is not 
easy to introduce for deducing the edge-like smearing 
function. Recent lunar profiles reconstructed from the 
altimetry data of the Kaguya mission (see Koutchmy 
et ah, 2013) can now be introduced, but it does not 
seem to have been used to deduce the PSF of the SOT 
of Hinode. Mathew et al. (2009) determined the PSF 
of the SOT also by analyzing the transit of Mercury. 
They introduced a combination of four Gaussians with 
different widths and weights in order to model the PSF. 
Gaussian functions can describe the narrow core of the 
PSF, however they cannot easily account for the far 
wings of the PSF describing the stray light (see for ex. 
Koutchmy & Koutchmy, 1974; Koutchmy, Koutchmy 
& Kotov, 1979 and the recent papers by Wedemeyer- 
Bohm 2008; DeForest, Martens & Wills-Davey, 2009; 
Wedemeyer-Bohm & Rouppe van der Voort, 2009), be¬ 
cause Gaussian functions decrease too fast with the ra¬ 
dial distance. In addition, matching the scattered light 
over a small range like the range given by the obser¬ 
vation of the transit of Mercury planet (10" diameter 
at time of the transit, which means that the wing of 
the PSF can be evaluated only for a short distance of 
order of the radius, e.g. 5") can produce only a limited 
knowledge of the PSF, especially beyond 5" radial dis¬ 
tance. The transit of the planet Venus of much larger 
diameter (60") was recently used in Yeo et al. (2014) 
for the analysis of the PSF and for deducing the stray 
light correction of the smaller in diameter SDO/HMI 
telescope. The stray light of the SOT is produced by 
the irradiation of the entrance of the telescope by the 
full Sun, although later the Gregorian mirror illumina¬ 
tion is limited to the field of view (FOV); ideally the 
far wings of the PSF should be determined, accounting 
for the value of the solar radius for components before 
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the Gregorian mirror and after, taking into account the 
FOV. 

Here we introduce a combination of a Gaussian (to 
represent the core of the PSF) and a Lorentzian func¬ 
tion (to represent the stray light at large distances) fol¬ 
lowing a suggestion by Johnson (1972). It includes the 
far wings of the PSF of the SOT, similarly to what was 
done in Koutchmy, Koutchmy and Kotov (1977), to 
deduce the corrected extreme limb darkening functions 
for the near infra-red wavelengths as measured with 
the vacuum telescope at Sacramento Peak (Richard B. 
Dunn Solar Telescope). Independently, a similar PSF 
was also proposed by DeForest, Martens, and Wills- 
Davey (2009) for describing the PSF of the space-borne 
TRACE telescope. We use our deduced PSF to decon¬ 
volve the blue continuum images of the SOT and will 
discuss several new results that were straightforwardly 
deduced. Removing scattered light considerably im¬ 
proves the RMS contrast of the granulation and this is 
critical for fine structures such as filigree and extreme- 
limb features and even more for studying bright dots 
and structures in sunspot umbrae and penumbrae, in¬ 
cluding their dynamical behavior. 

2 Observations 

We first used images taken during the Venus transit on 
6th July 2012 recorded with the broadband filter imager 
(BFI) of the solar optical telescope (SOT) onboard the 
Hinode satellite (Tsuneta et al, 2008; Suematsu et ah, 
2008) to determine the Point Spread Function (PSF) of 
the telescope, (see Fig. 1 for a selected image). For the 
sake of simplicity in interpreting the photometric data, 
we restrict our analysis by avoiding images that were 
taken with filters containing contributions from chro¬ 
mospheric lines or from a molecular band that need 
a special treatment. The blue continuum images are 
selected because they correspond to the shortest wave¬ 
length available, presumably giving the best spatial res¬ 
olution because of the favorable ratio with the given 
value of the aperture of the telescope. In addition, 
this filter gives the best temperature contrast among 
all continuum filters available. Similar approach was 
recently used by Louis et al. 2012 in their studies of 
umbral dots. The filter has a central wavelength at 
4504.5 A and a band width of 4 A corresponding to a 
low line-blanketing effect. It is an excellent window of 
the continuum solar spectrum over the whole solar sur¬ 
face, see for ex. the recent solar atlas of Fathivavsari, 
Adjabshirizadeh and Koutchmy (2014). 

Well focused images were selected, and each pixel 
corresponds to 0.05448 arcsec on the surface of the Sun. 


Exposure time is 205 ms, the cadence is 60 seconds and 
the field of view is approximately 111x111 arcsec ^.We 
also used a set of images taken on 1th March 2007 in 
blue continuum filter of a single unipolar almost sym¬ 
metric sunspot images to make our movie. 125 images 
taken between 00:14:51 and 00:59:53 were selected in or¬ 
der to make it; the exposure time equals to 102 ms and 
the cadence varies between 6 and 60 seconds. Each pixel 
corresponds to 0".05448 and the field of view (EOV) is 
approximately 56x28 arcsec Eurthermore, an image 
of this single sunspot viewed close to the disc center but 
with a larger EOV (223x111 arcseconds) is also used. 
It is observed one day before on 28th Eebruary 2007 at 
the same wavelength and with the same pixel value (ex¬ 
posure time equals to 77 ms) to compute RMS contrast 
at disc center for magnetic and nonmagnetic regions. 

Einally, we qualitatively use a Gall H line filter im¬ 
age (Exp time 0.205 ms with the same pixel value, 
0."05448) taken at the same time in order to charac¬ 
terize these regions in details. All images were first 
corrected for dark current, fiat field and bias using the 
fg_prep routine from the Solar Software. 

3 Data processing and deduction of the PSF 

The 10-arcsec-diameter Mercury disc was already used 
to compute and model the PSE (Mathew, Zakharov, 
Solanki, 2009; Wedemeyer-Bohm, 2008). The Venus 
disk gives a much better evaluation for the wings of the 
PSE because its radius is approximately 6 times bigger 
than the Mercury radius during the transit across the 
disk of the Sun. Scattered light can be seen inside the 
Venus disc up to a larger distance in comparison with 
Mercury disc (Yeo et ah, 2013). In order to make the 
radial profile of Venus, after making several tests of 
the resulting signal/noise ratio (S/N), we decided to 
use the averaging of 30 rows (15 rows over and 15 rows 
under the apparent center of the planet) with respect to 
the diameter without significantly losing the resolution 
while the small curvature of the edge of Venus keeps 
the edge within half a pixel. Eig. 1 shows an observed 
image during the Venus transit and Eig. 2 shows the 
average radial profile of Venus using 30 rows. Because 
of the limited size of Venus, we further decided to use 
the limb of the Sun and the available part of the image 
of the sky to compute the far wings of the PSE taking 
into account the light recorded from far distance outside 
the disk. 

In order to make a radial profile of the solar limb for 
deducing the smearing function, we improve the signal 
to noise ratio by averaging many rows, taking into ac¬ 
count the curvature of the limb as was already done 
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Fig. 1 Observed SOT image during the Venus transit using 
the blue filter. Note the limb darkening on the partial image 
of the Sun. 


in the case of the profile of the Venus planet. Now, 
the curvature of the limb is much smaller and also we 
can go much further out from the sun where the sig¬ 
nal is weak and the instrumental aureole is observed 
as in Koutchmy, Koutchmy and Kotov (1974). We can 
average up to approximately 400 rows without the cur¬ 
vature of the limb producing an effect surpassing half 
of one pixel. Indeed making derivative profiles of the 
solar limb for different numbers of rows (in pixel) par¬ 
allel to a selected diameter direction (200, 400 and 600 
rows), the resulting derivative profiles of the solar limb 
are evaluated in terms of ” smoothing” when more pixels 
are added but taking into account the effect of curva¬ 
ture. It can be seen that when averaging a high number 
(above 400) the derivative of the limb is broadened. We 
use an averaging over 200 rows that corresponds to an 
effect due to the curvature of the solar limb of order of 
0.14 pixel (half value of the sagitta). Fig. 3 shows the 
average limb profile of the Sun for a 200 row averaging. 
The ripples appearing at the left hand side are mainly 
due to the disk granules and heterogeneities at extreme 
limb (we will discuss this effect further). 

Fig. 4 shows the average derivative profile of the so¬ 
lar limb for 200 rows. The left part of Fig. 4 is rough 
because of the influence of the granules but the right 
side does not show the influence of granules, so this 
curve should better be used. We then cut it through its 
center of gravity to make it symmetrical. We can fit the 
resulted curve with an appropriate function which is the 
Abel Transform (AT) of the PSF, assuming the extreme 
limb is a step function at least at the scale we consider 
here, based on the theoretical expectation we have of 
the solar limb from a hydrostatic ID model (Vernazza, 
Avrett and Loeser, 1976). In other words, the limb 
derivative after the inflection point is the projection 
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Fig. 2 Average radial profile of the limb of Venus using 
a 30 rows average from Fig. 1. Inserted is the cross sec¬ 
tion of the Mercury planet limb from Mathew, Zakharov, 
Solanki, 2009. Note the difference of scales in arcsec with 
the comparable level of stray light over the shadow of the 
planets. 


of a circularly symmetric PSF along a set of parallel 
lines and represents the Line spread function (LSF), 
(see Aime, 2007 for a more elaborate discussion). The 
PSF could be deduced from the LSF by using an inverse 
Abel transform, see Bracewell (1976). After finding the 
center of gravity (which is very close to abscise of the 
maximum intensity) we cut the peak from this point 
and reverse the right side to replace with the left side 
and have a symmetrical curve see Fig. 5. Before fit¬ 
ting the resulted curve, we found that it is better to 
rebin the data by a factor greater than one (for exam¬ 
ple 3 or 4) to reduce the role of noise during the fitting 
procedure and also to go further out, especially in far 
wings. We now found that the PSF can be approxi¬ 
mated using the sum of a Lorentzian plus a Gaussian 
function. It means that the Lorentzian function repre¬ 
sents the stray light that comes from far distances while 
the Gaussian function corresponds to short distances, 
including the diffraction by the entrance aperture, close 
to the Airy function. The Abel transform being addi¬ 
tive, we can fit this curve with the linear combination 
of the Abel transform of the Lorentzian function plus a 
Gaussian function. The Abel transform arises when cir¬ 
cularly symmetric distributions in two dimensions are 
integrated along one dimension. The Abel transform 
A[f(x)] of the function f(r) is defined as 


r-\-oo 

[A[fix)] = 2 / 

J X 


f{r)r dr 
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Fig. 3 Radial profile of the limb of the sun using a 200 
rows averaging. 
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Fig. 5 Full symmetric derivative profile of the solar limb 
after rebinning the data by a factor 4 to show it further out. 
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Fig. 4 Derivative of the radial profile of the limb of the 
Sun shown in Fig. 3. 
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Fig. 6 Symmetric radial profile of the solar limb after re¬ 
binning by a factor 4 (black) and the fitted line (red) 


This transform is additive (Bracewell, 1976), it 
means that 


^[f{x) + g{x)] = A[f{x)] + A[g{x)] (2) 

To derive the PSF, we fit Fig. 5 with the Abel trans¬ 
form of the Lorentzian plus the Gaussian functions us¬ 
ing weights pi; we call it R: 

R = A{poL + psG) = poA{G) + psAiL) (3) 


where L and G are Lorentzian and Gaussian function, 
respectively. This can be written as: 


R = 


{p\ 


Po'X _ 

{X-P2)^)^ 


-(x-p2)^ 

+ P3V^P4e ^’’4 


(4) 


where po and ps are coefficients (weight) that indi¬ 
cate the relative contributions of the Gaussian and of 
the Lorentzian functions in the fitting procedure, p 2 is 
the location of maximum (Fig. 6). The result of the fit 
is depicted in Fig. 6, and its parameter are as follows: 


Po = 3.144, Pi = 1.395, p2 = 0, P3 = 1.088, p4 = 2.8 (5) 


As discussed before, equation (1) is the Abel trans¬ 
form of the PSF which is a linear combination of the 
Lorentzian and the Gaussian functions. 

Un -(x-P2)^ 

PSF{x) = , ^ + pse (6) 

Pi - {x-p2y 

Goefficients po lo p 4 are obtained from the fit to the 
observed functions. In Fig. 7 and Fig. 8 three dimen- 
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Fig. 7 Three dimensional plot of our obtained PSF, in 
logarithmic scale 

sional plots of our PSF is shown; the Mathew et al. 
(2009) result is also shown to make a qualitative com¬ 
parison. Fig. 9, and Fig. 10 show them together for a 
better quantitative evaluation. 

As it can easily be seen from Fig. 9, that the Mathew 
et al. (2009) PSF goes rapidly to low values after 3”, 
so it could not fully take into account the stray light 
coming from moderate and great distances (when con¬ 
volving using the PSF, a 920” radius Sun should in 
principle be considered, depending on the direction and 
on the position of the feature at the solar surface). Fur¬ 
thermore it does not have a physically and meaningful 
shape, since below 10“5 (or 6” in radial distance) its 
slope changes suddenly instead of decreasing slowly and 
continuously. They did not elaborate on the reason for 
this sudden change in the slope. As far as the core of 
the PSF is concerned, our PSF is rather narrower in the 
top part (Fig. 10) and this is what is expected for an 
aperture with a central obstruction, as was noticed from 
the modulation transfer function concept discussed by 
Wedemeyer-Bdhm and Rouppe van der Voort (2009). 
Accordingly, we believe the PSF we propose is a good 
approximation of the effective PSF of the SOT, at least 
in the frame of the assumption of an isotropic and shift- 
invariant PSF. Fig. 9 and 10 show also the Airy pat¬ 
tern for a circular homogenously illuminated aperture 
of diameter d=50 cm without a central obstruction. To 
prove the effectiveness of our PSF in removing the stray 
light, we show in Fig. 11 and Fig. 13 some radial pro¬ 
file near the limb of the Sun, taken at a position outside 
of any sign of activity. Only fluctuations at the gran- 


Fig. 8 Three dimensional plot of Mathew, Zakharov, 
Solanki, 2009 PSF in logarithmic scale. 

ulation scale are seen. Noticeably, no light is observed 
outside the limb after removing the effect of the PSF 
(after deconvolution, see Fig. 11 and Fig. 15) and the 
inflection point usually considered as showing the limit 
of the solar diameter, after corrections, is slightly inside 
the observed inflection point. 

4 Applications to different solar photospheric 
phenomena 

Restoration of images was readily performed using the 
method already successfully used for SOT data by Louis 
et al, 2012. It takes advantage of the available IDL 
Max-likelihood routine (Richardson, 1972; Lucy, 1974) 
from the Astrolib which deconvolves the observed im¬ 
ages provided the instrumental PSF is known or as¬ 
sumed. This routine performs iterations based on the 
Maximum Likelihood solution for the restoration of a 
blurred image. The routine does iterations and updates 
the current estimate of the image by the product of the 
previous deconvolution and the correlation between re¬ 
convolution of the subsequent image and the instrument 
PSF (Richardson, 1972; Lucy, 1974; Mathew, Zakharov 
and Solanki, 2009; Louis et ah, 2012). A parameter that 
we introduce when using this routine is the number of 
iterations to use, which depends on the quality of the 
original data. We used an empirical method instead of 
trying to discuss the relevance of the spectral distribu¬ 
tion of the noise present in the data for optimizing the 
number of iterations to be used. Indeed, we are not 
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Fig. 9 In log scale i) the theoretical Aity pattern for a full 
50 cm aperture and a 450 nm wavelength (in grey), ii) the 
deduced effective PSF used in this paper (black full line) 
and iii) the one used by Mathew, Zakharov, Solanki, 2009 
(dashed line) 


Fig. 10 In linear scale, the Airy pattern for a fully illu¬ 
minated 50 cm aperture at 450.4 nm (gray solid line), our 
deduced effective PSF (black solid line) and the PSF used 
by Mathew, Zakharov, Solanki, 2009 (dashed line). 


really limited in using a large number of iterations in 
order to increase the contrast and, accordingly, the vis¬ 
ibility of the features. We do not think that more than 
several hundreds iterations would help, as suggested by 
the results of the following test. 

A granulation field situated near the single sunspot 
we studied (see further), was selected to be representa¬ 
tive of both the quiet Sun and of the magnetic granu¬ 
lation inserted inside (see Baudin, Molowny-Horas and 
Koutchmy, 1997 for a description of the magnetic (MG) 
and the non-magnetic granulation (NMG) that we will 
analyze further out). We used this test field to look 
how the deconvolution works by making 50, 100, 200, 
300, 500, 600 ,700 and 800 iterations; Fig. 12 shows 
the result, to be compared with the original image in¬ 
serted at the center. The impression is that the more 
iterations you use, the more small scale structure we 
see, the obvious limitation being the structures of non 
solar origin artificially produced by the noise, with 50 
iterations we already see all features of the intergranu¬ 
lar lanes. We could decide to stop the iteration process 
when spurious structures appear based on a qualitative 
reproducibility criteria when using consecutive in time 
frames. A similar procedure was used in the Koutchmy 
(1977) paper dealing with the deconvolution of a filigree 
image using the Weiner filter. This is especially true 
in the part where some magnetic field was measured 
with corresponding very fine structure called abnormal 
granulation (Dunn and Zirker, 1973), as at the bottom 
right of each box of Fig. 12. This abnormal granula¬ 
tion is also now called the magnetic granulation (MG), 
e.g. Baudin, Molowny- Boras and Koutchmy (1997). 
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Fig. 11 Observed (full line in black for the averaged pro¬ 
file) and restored limb (grey lines for narrow scans and in full 
red line for the averaged scan) of the Sun, vertical dashed 
line shows cos^ =0.1. 

Without going too far into the details, we note that 
recent numerical simulations of the granulation are ca¬ 
pable of reproducing many details of the granulation 
(Wedemeyer-Bdhm and Rouppe van der Voort, 2009). 

In order to make a more quantitative evaluation of 
our results, we computed the RMS of intensity fluc¬ 
tuations for each deconvolved image and finally, the 
chi-square value of the successive differences between 
the resulting image re-convolved with our PSF and the 
successive deconvolved images with different numbers 
of iterations: 
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Fig. 12 Observed (Mosaic made of the same granulation 
field after deconvolving using different number of iterations 
(see the number put at top left in each case). The original 
image is at the center and it is covering 25^^x25^^area. 


Fig. 13 Deconvolved limb of the Sun in different positions 
across the limb (gray lines) as shown in Fig. 14, and their 
average (red and blue lines) for observations taken 6 minutes 
apart. 


iterationnumber 

RMS Contrast 

Chisquare 

12 

0.213 

317.31 

25 

0.236 

302.43 

50 

0.250 

294.36 

100 

0.261 

286.06 

200 

0.273 

283.88 

300 

0.282 

283.82 

400 

0.289 

284.02 

500 

0.295 

284.03 

600 

0.300 

283.8 

700 

0.305 

283.40 

800 

0.309 

282.93 

Original 

0.109 



From the examination of the chi-square values, it 
seems that by using 200 iterations, we reach a rela¬ 
tive minimum in the and further, some fluctuation 
is seen. Accordingly, this is the maximum number of 
iterations adopted for the deconvolution of our images, 
including the case of the sunspot movie we consider be¬ 
low. As far as the value of the RMS of the granulation 
is concerned, with 50 iteration we already reach a value 
in agreement with Wedemeyer-Bdhm and Rouppe van 
der Voort (2009) and the synthetic values reported by 
the authors. We recall that the deconvolution is done 
over a limited FOV; we do not take into account the 
whole Sun for removing the stray light which is partly 
reduced by the Gregory optical system used for the SOT 
observations. Indeed a more extended analysis could be 
done on this topic, but it is outside the scope of this pa¬ 
per, see Wedemeyer-Bdhm and Rouppe van der Voort 
(2009) for a discussion of results coming from deeply de¬ 
convolved images when taking into account the noise. 


4.1 Study of the extreme limb 

ID models of the solar atmosphere can be tested by 
the behavior of their predicted center-limb variation of 
intensities in the continuum. The test is particularly in¬ 
teresting at the extreme-limb where more outer layers 
are concerned because line of sight integration effect. 
Then, the influence of the emerging magnetic field can 
be important, at least where the magnetic flux tubes fan 
out, as in the network, and are assumed to be respon¬ 
sible for the heating of the chromosphere. Recall that 
the visible limb of the Sun is typically situated at a dis¬ 
tance of the order of 370 km from the reference height in 
photospheric models T 5 = 1 , which means in the vicin¬ 
ity of the temperature minimum region, below or above 
depending on the model atmosphere (see Rutten, 2002 
for an extensive review of historical and recent models 
of the Solar Atmosphere). Additionally, the question 
of the solar diameter variation and its shape, and the 
extreme-limb of the Sun in its most relevant detail re¬ 
cently received a lot of attention (see for ex. Rozelot 
and Fazel 2013), not to mention the linear polarization 
properties at the extreme-limb which is another way to 
study this region. The topic is treated in so many pa¬ 
pers including the so-called temperature bi-function ef¬ 
fect (Rutten, 2002) that we restrict our discussion here 
to the analysis of measurements at 450 nm and at dis¬ 
tances from the classical limb, defined by the inflection 
point, well under 1 ”. 
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Fig. 14 Selected image of the limb of the Sun where 
extreme- limb scans of Fig. 13 were obtained. The boxes 
show the extension of each individual scan used to get an 
average; the positions of the boxes correspond to the range 
covered by different successive averages, shown by grey lines 
in Fig. 13 by shifting the boxes along the limb. 

We first examine in Fig. 13 the averaged radial pro¬ 
files of the limb of the Sun after deconvolution. Fluc¬ 
tuations are clearly apparent, up to the extreme limb, 
including a very small limb brightening, whose ampli¬ 
tude diminishes with the averaging procedure (see Fig. 
14 the range of the limb positions we used). A more ex¬ 
tended analysis is probably needed to give the precise 
value which possibly shows a pole equator variation. 
To further investigate the extreme limb, we took the 
derivative of the deconvolved averaged variations, see 
Fig. 15 and compared it with the theoretical results 
available from the literature, see Fig. 16. We found an 
excellent agreement, especially on the ’’width” of the 
limb, although inside the limb, our small extreme limb 
brightening is not reproduced in the former theoreti¬ 
cal analysis of Lites, 1982 of the solar limb for different 
continuum windows, based on the use of the VAL classi¬ 
cal model (Vernazza, Avrett and Loeser, 1976) that did 
not yet consider the so-called ’’bifurcation” of temper¬ 
ature models (see Rutten, 2002) . The FWHM of the 
derivative of the theoretical limb (0.217 arcsec) is in ex¬ 
cellent agreement with the FWHM of our deconvolved 
averaged limb profile (0.214 arcsec). Note that we took 
the Lites, (1982) data after interpolating his results to 
get the behavior exactly at 450 nm. Some more recent 
values of the limb darkening intensities were published 
in the literature in the frame of discussions of the solar 
diameter determinations with different methods, but 
usually they do not reach the spatial resolution of the 
Lites (1983) work such that we do not mention them. 
We recall that the Lites (1983) analysis was done to 


18 19 20 21 22 



Fig. 15 Extreme solar limb averaged variations (greatly 
magnified) after deconvolution in relative units, similar to 
the units of Fig. 15. Note that the stray light outside the 
limb is completely removed. 

look at the so called ’’rugosity” of the limb, which are 
small-scale variations of amplitude in the height of the 
solar surface, presumably due to a solar granulation 
thought to be present in the higher layers where today 
we know that some reversing of the contrast occurs. 
The classical granulation recorded over the solar sur¬ 
face with high contrast could show a signature near the 
limb. The contrast is supposed to reverse at a height of 
the order of 150-200 km above rs = 1 from models (see 
for ex. Gadun, 1995 and observationally, Koutchmy, 
1994), at least from the analysis of images taken simul¬ 
taneously in the continuum and in photospheric lines 
formed near these heights. It means the reversing oc¬ 
curs at distance from the limb where T 5 is lower than 
0.1 and the effective temperature less than 5000K, from 
the VAL model. Numerical 3D simulations of the so¬ 
lar granulation published for the last decade beautifully 
confirm this 1st approximation behavior. Fig. 17 is a 
greatly magnified part of the solar limb taken from our 
deconvolved image used to build Fig. 13-16 to show the 
fluctuations. 

Finally we note that the question of the extreme limb 
brightening is an old challenge for observers. The best 
method for observing this effect is attempted at solar 
total eclipses (Bazin and Koutchmy, 2013) because the 
level of stray light is then drastically reduced. Rosen 
and Poss (1982) reported a limb brightening in the last 
1 arcsec from the infiection point taken as the position 
of the limb). We do not find an agreement with our 
value that shows the limb brightening much closer to 
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Fig. 16 Theoretical extreme limb variation (dashed line) 
from Lites, 1982 and its derivative (cross symbol). Solid line 
shows the derivative of the deconvolved limb of Fig. 14. 



Fig. 17 Greatly magnihed image of the extreme limb vari¬ 
ations in negative (the sky outside the bright solar disk is 
white). Note the structure observed up to the extreme limb 
and the very narrow and faint irregular limb brightening. 

the limb (see Fig. 13 and qualitatively using the greatly 
magnified image of the limb shown in Fig. 17.) Eclipse 
data based on the analysis of a recording of the eclipse 
light curve probably suffers from the very irregular edge 
of the crescent of the Moon that is difficult to take 
into account. Using fast flash spectra or fast slit-less 
spectra, a much better analysis is possible, because the 
edge of the Moon is imaged, see Bazin and Koutchmy 
(2013), from their reported flash spectra, we noticed 
that the rather strong ’Tow FIP” emission line of Till is 
close (450.18 nm) to the center wavelength of the blue 
continuum filter used in the SOT observations. This 
line could contribute at the extreme limb and possibly 
explain the observed small limb brightening. 


4.2 Analysis of the Granulation contrast over the 
surface of the Sun 

We applied the deconvolution procedure using our PSF 
on the granulation observed with the excellent focus 
of images taken with a pixel size of 0.05448 arcsecond 
during the transit of Venus. We evidently applied this 
deconvolution procedure in order to restore the image 
of Venus during its transit near the egress of the planet. 
The result gives us confidence in our applied correction, 
as all the spurious light observed on the dark disk of the 
planet is removed but the study of the aureole of Venus 
is clearly outside the scope of this paper and will be 
considered later. 

The deconvolution applied to the granulation of the 
quiet and the less quiet Sun immediately brought inter¬ 
esting results. The normalized RMS contrast in quiet 
photosphere area drastically increases after the correc¬ 
tion is applied. The normalized RMS contrast at disk 
center in quiet photosphere is typically 0.11 before cor¬ 
rection (using original observations) and it is 0.25, using 
50 iterations, see Fig. 12 and further. The agreement 
is excellent compared to the results reported in the ex¬ 
tended analysis of Wedemeyer-Bohm and Rouppe van 
der Voort, L. (2009) who developed a different method 
of correction, based on assumptions quite similar to 
ours and a PSF described with another function. It 
makes these results regarding the high values of the 
granulation contrast after correction for the PSF, rather 
solid and they are now a valuable test of theoretical 
models of granulation. 

We will not extend our analysis of the granulation 
too far because we basically agree with the results al¬ 
ready exposed in the exhaustive Wedemeyer-Bohm and 
Rouppe van der Voort (2009) paper. We prefer to in¬ 
vite the interested reader to study that paper that in¬ 
cludes a theoretical part taking advantage of synthetic 
images produced by these authors from up to date com¬ 
puter simulations and models. We just add here in 
short some quantitative results that seem to us still in¬ 
teresting to show. It is in particular about the near 
limb granulation. Note that the study is made with¬ 
out taking into account the influence of p-mode inten¬ 
sity oscillations that indeed near the limb should not 
significantly change the results, from the discussion in 
Baudin, Molowny- Horas and Koutchmy (1997). Fig. 
18 shows our normalized RMS contrast before and af¬ 
ter deconvolution, for different distances from the limb 
(inflection point); typical values of the view angle 0 
from the normal to the local surface defined as usual 
with cosO = /i are also shown. As expected, the RMS 
decreases towards the limb, similarly to the limb dark¬ 
ening. The decrease is partly due to the foreshortening 








11 




Fig. 18 Variations of the RMS contrast near the solar 
limb, from both the original image and from the corrected 
image reflecting the behavior of the granulation contrast. 
Vertical dashed lines show cos 0 =0.1, 0.2, 0.3 respectively 
from left to right. 

effect which makes narrower granules in the direction 
towards the limb. Due to the effect of dispersion of the 
granule intensity variations, the RMS contrast reflects 
the limitation of the amount of data used, however, 
no clear relative brightening when going to the limb is 
observed, noting that the part of the image we used did 
not show any evidence of faculae, see Fig. 17. 

Indeed, an interesting feature of granulation is the 
influence of the very small scale magnetic elements orig¬ 
inally called filigrees that are more important in plages 
and also going to the limb. This magnetic granulation 
(MG) was carefully compared to the ’’non-magnetic” 
granulation (NMG) by Baudin, Molowny- Horas and 
Koutchmy 1997. Unfortunately, they used rather small 
spatial resolution and not corrected ground-based ob¬ 
servations, leading to possibly limited results. Using 
a larger field image from Febr. 28, 2007 coming from 
the observation of the same single sunspot region that 
we discuss after in 4.3, we produced a corrected image 
using our PSF. We then computed the RMS of granu¬ 
lation for both a non- magnetic region granulation field 
NMG and for a magnetic one MG, see Fig. 20. We 
found, after deconvolution using 200 iteration(a bigger 
number of iteration were used in order to increase the 
contrast and the visibility of small scale features): RMS 
(MR) = 0.297; RMS (NMG) = 0.295. The correspond¬ 
ing values for the original image (not corrected) RMS 
are 0.123 and 0.122. Finally, we show in Fig. 19 the 
results of making a histogram analysis in this region 
taken near the solar center, see Fig. 20, to illustrate 
the skewness of the brightness distribution, similarly to 


Fig. 19 Histograms of the distribution of brightness ob¬ 
served over the deconvolved image of the granulation. In 
red, for the magnetic granulation; in black for the nonmag¬ 
netic granulation. 

what was discussed in Koutchmy, 1994 from the anal¬ 
ysis of the IR (e.g. deeper) granulation observed in 
the opacity minimum region. The agreement is excel¬ 
lent with what was obtained by Wedemeyer-Bohm and 
Rouppe van der Voort, L. (2009) for the deconvolved 
normal granulation image (see their Fig.5, left part), 
taking into account their scale that uses the logarith¬ 
mic of intensities. Making an evaluation of the differ¬ 
ence corresponding to this bi-modal distribution which 
presumably corresponds to a set of bright and to a set 
of dark features that can be translated in temperature, 
we see an overall agreement with the IR results that 
gave a ratio 1: 0.7 of dark areas surpassing the bright 
areas. However now, the difference in temperature is 
larger because IR data were not corrected for smearing. 
In addition, this discussion does not include the effect 
of filigrees predominantly inserted in the dark lanes of 
the MG. It needs in a finer analysis. Fig. 20 clearly 
shows the difference of the field of MG and the field 
of NMG selected using a SOT image taken simultane¬ 
ously with the H line Gall filter where plage filigrees 
are considerably enhanced. Fig. 19 compares what is 
seen in the quiet part of the image (NMG) with what 
is seen in a magnetic granulation (MG) region using 
the histogram analysis. The histograms look different 
around the average values of intensities (central part of 
the histogram); a careful examination shows a definite 
but really small global shift of the MG histogram to¬ 
wards higher values of intensities. We finally note that 
the skewness of the histograms revealing a bi-modal 
distribution, is more clearly seen for the NMG. 
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Fig. 20 Image taken in the H line of Call on Febr. 28, 
2007 showing the selected regions for comparing the mag¬ 
netic granulation MG (box at left) and the nonmagnetic 
granulation NMG, box at right bottom part. 

4.3 Analysis of a single sunspot 

We now report some preliminary results correspond¬ 
ing to the main motivation of our study, the analysis 
of sunspots with the aim of understanding their origin 
and evolution (Solov’ev and Kirichek, 2014). Correc¬ 
tions needed to apply to the data play an even more 
important role than in the case of the quiet solar sur¬ 
face features, because the true umbral intensity of a 
sunspot core is very low. Accordingly, the stray light 
level in the umbra is frequently at a comparable level to 
the umbral brightness (see for ex. Adjabshirzadeh and 
Koutchmy, 1983). Moreover, the umbrae of sunspots 
show tiny bright features called umbral dots (UDs), of 
very small diameter. In the early 80s it was shown 
that their size is close to the diffraction limit of a 76 
cm diameter telescope in the blue region (see Adjab¬ 
shirzadeh, and Koutchmy, 1982-83), which means typ¬ 
ically 0.”25 (180 km), with a distribution over the um¬ 
bra that makes them difficult to be definitely charac¬ 
terized. Indeed, clumps of UDs are seen, as well as 
alignments of UDs into elongated lanes. UDs are dy¬ 
namical objects and their proper motion was studied 
(see Molowny-Horas, 1994). Qualitatively, proper mo¬ 
tions can be demonstrated using a movie. This dynami¬ 
cal behavior and time variations are difficult to analyze 
from ground-based observations where the quality of 
the seeing is variable; it is much more convenient to 
analyze sequences with seeing-free high resolution cor¬ 
rected images. 

Louis et al. (2012) published an extended analysis of 
UDs from SOT blue filtergrams processed in a way sim¬ 
ilar to our observations. However they used corrections 
with the PSF of Mathew, Zakharov and Solanki, (2009) 
which is quite different from the one we deduced here 


(Fig. 9 and Fig. 10), so that we anticipate differences 
in the results. The Louis et al. (2012) paper already 
describes in details the observations of March 1st 2007 
that we used. Here we first provide some rather quali¬ 
tative and preliminary results that apparently have not 
previously been seen, to our knowledge. 

To evaluate the effect of the restoration method 
with our PSF, we apply it to the deconvolution of this 
’’known” sunspot of March 1st (2007). Fig. 21 (top) 
shows an observed image and Fig. 21 (bottom) shows 
its deconvolved image. It can be seen from the exam¬ 
ination of these images that the deconvolution ’’brings 
out” some umbral fine structure , such as UDs, with¬ 
out the need of over-exposing the picture, which means 
that they reach intensity comparable to the intensity of 
the surrounding photosphere. However, we note that 
the dots are observed to have a wide range of intensities 
and the corrected contrast depends on the number of it¬ 
eration used. Regarding the peripheral UDs, their cor¬ 
rected intensities are even more impressive (as it was al¬ 
ready reported in the Sobotka and Puschmann, (2009) 
images from the 1-m Swedish Solar Telescope reaching 
a 0”14 resolution). We finally decided to use 200 iter¬ 
ation for improving the visibility of all umbral features 
noting that the deconvolution does not produce any ar¬ 
tifact effect like the undesirable ringing effect around 
the smallest point-like bright dots. 

The contrast of penumbral filaments and the mag¬ 
netic granulation around sunspot are also greatly im¬ 
proved. A movie was prepared from deconvolved blue 
continuum images taken between 00:14 and 00:59 (125 
frames in 45 minutes) to see the dynamical behavior 
of fine structures inside umbra and penumbra. This 
is readily possible because of the absence of seeing 
in space-borne observations. Regarding the dynam¬ 
ics of the penumbra, the movie made of corrected 
frames gives some qualitative evidence of rather large 
scale emerging penumbral diffuse dark radial magnetic 
strands moving together with elongated bright quasi- 
radial filaments with evidence of twist and dark lanes 
inside. Simultaneously, there are bright features mov¬ 
ing inwards with detachments making the peripheral 
UDs. At smaller scale, a more careful analysis, possibly 
with additional spectroscopic material, should be done 
before discussing the physics of the penumbra. The 
movie however shows excellent examples of the inwards 
motion of bright rather elongated elements towards the 
periphery of the umbra produces peripheral UDs with a 
definite excess of brightening, suggesting that nonlinear 
processes producing a shock could occur at the umbral 
boundary where the magnetic field sharply changes of 
direction. In addition, the penetration process into the 
umbra continues with evidence of splitting. These phe- 
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Fig. 21 Observed SOT image of the single sunspot of 2007 
March 1 (top) and the deconvolved image 

nomena will be addressed more quantitatively in the 
future. 

A last result we want to show concerns the most crit¬ 
ical part of the sunspot, from the point of view of the 
origin of the strong magnetic flux emerging from the 
deep layers responsible for a large part of the activity 
of the Sun which is concentrated in the core. Fig. 22 is 
a display made to show the structure of the core, with 
UDs visible inside, including i) an alignment pattern of 
UDs partly making a very weak bridge, with a dark lane 
inside, not far from the central part of the core; their 
morphology reminds us of what is observed in the clas¬ 
sical light bridges of sunspots (Secchi, 1870) but here, 
we deal with a very low intensity feature; ii) an area 
in the core close to the weak bridge, where, surpris¬ 
ingly, no UDs are present, making a sort of flat area of 
approximately 3"x2" surrounded by several elongated 
UDs, corresponding to the darkest part of the umbra. 
We have not found in the literature a description of a 
similar feature in the umbrae of sunspots. These pecu¬ 
liarities cannot be detected on the uncorrected images. 
We believe that further analysis is needed to confirm 
their presence before speculating on the physical signif¬ 
icance of this finding. Hopefully, such confirmation will 
come from a much deeper analysis of our movie ma¬ 
terial permitting to tackle very small scale dynamical 
phenomena. 



Fig. 22 Negative display in both linear and logarithmic 
scales of intensity modulations of the core of the sunspot 
showing peculiar features of the umbra. 

5 Conclusion 

The blue continuum high resolution images of the SOT 
of the Hinode mission were corrected for stray light 
based on several observations, including the use of im¬ 
ages of the transit of planet Venus and the case of the 
extreme limb ID variation. A combination of Gaus¬ 
sian and Lorentzian Functions were used to construct a 
new telescope PSF. This combination of functions well 
represents the diffracted light that comes from narrow 
angles (core of PSF) and the stray light that comes from 
large angles (far wings of the PSF). 

The subsequent deconvolution using iterations of im¬ 
ages greatly improves the contrast of any image of the 
solar surface. In the case of the aureole of the planet 
Venus, an impressive gain is obtained for looking at 
its real intensity variations, it removes all the spurious 
light from the dark disk of the planet and it gave us 
confidence in our method of restoration. 

The normalized RMS contrast of the quiet photo¬ 
sphere is increased by more than a factor 2 and some 
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statistical properties of the granulation, inside and out¬ 
side magnetic regions such as a plage with filigrees in¬ 
serted, are deduced. Furthermore the corrected images 
show a new result regarding the extreme limb of the 
Sun, including a faint limb brightening at 0”2 of the 
edge of the Sun where indeed elongated along the limb 
irregularities are still observed. No interpretation is of¬ 
fered but it suggests that the simultaneous analysis of 
images obtained with different filters could help. 

Finally, deconvolved images of a sunspot also show 
some new dynamical penumbral structures and espe¬ 
cially some new umbral structures much more clearly, 
such as aligned umbral dots (UDs) seemingly making 
the start of an umbral bridge. The core of the um¬ 
bra also shows an area where UDs are absent, which 
seems difficult to reconcile with an assumed vigor¬ 
ous magneto-convection present permanently inside the 
umbra. The deduced movie is proposed as additional 
material for the future analysis of small scale dynam¬ 
ical phenomena. It needs careful study, including cor¬ 
relation analysis of successive images with short time 
interval that we intend to describe in a forthcoming pa¬ 
per. We believe this is now possible because i) images 
are free of any turbulent effect, ii) their photometric 
properties are excellent and hi) the proposed correc¬ 
tion based on the use of a rather simple PSF makes 
the deduced results solid, confirming that we deal with 
wonderful material. 
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